The CNO Bi-Cycle 


A second way of processing 'H into ^He is through a series of nuclear 
reactions involving the different species of Carbon, Nitrogen, and 
Oxygen. The principle reactions are as follows: 


120 +H — BN +y (14.1.1) 
BN — DOC et + v, (14.1.2) 
130 +H — MN +y (14.1.3) 
MN + 1H — 150 +y (14.1.4) 
150 — IN + et + re (14.1.5) 
15N + 1H — C+ *He (14.1.6) 


or 
15N +H — O +y (14.1.7 
160 + 1H — "F +y (14.1.8 
WEF — O +e +re (14.1.9 


170 + 1H — MN + *He (14.1.10 
Note that the net sum of these reactions is 
DO du TIE a eC ee 28 dep 


In other words, four protons are converted into two positrons, two 
electron neutrinos, and one “He nucleus; !?C serves only as a catalyst. 
Moreover, since the cycle can begin with any of the intermediate 
produces, it will occur if any of them exists. 





The differential equations for the full CNO bi-cyle are below. To sim- 
plify the notation, the number density symbol, N has been omitted. 
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where 7 is the lifetime of the given species against reactions with a 
proton (r = 1/AfNig), and Tg is the species lifetime against beta 
(or positron) decay. In the CNO cycle, all the reactions, with the 
exception of !"O(p, ^He)4N, are believed to be non-resonant. 








There are several items of note in the previous equations. First, 
when !?N combines with a proton, the result is an excited compound 
nucleus, !éO". This nucleus then decays into either normal !6O or 
into ^C + ^He. The probability of the former reaction depends on the 
initial energy of !9O, which, of course, depends on the energy of the 
reaction. Thus, the branching ratio, y, which appears in (14.1.11) 
and (14.1.17), is actually a mean value that is averaged over the 
Maxwellian velocity distribution. It is thus dependent (weakly) on 
temperature. Its value, y ~ 4 x 10-7, indicates that most CNO 
processing occurs through the first upper part of the bi-cycle. 


A second item to note is that !?N, 150, and !*F can be lost either 
by beta decay or by reactions with a proton. In all cases, however, 
the timescales for beta decay are much shorter than that for nuclear 
fusion (typically seconds vs. years). Thus, unless the temperatures 
are extremely high (T > 108 K), the fusion terms in these equations 
can be dropped. 


[Note: in the late stages of stellar evolution, temperatures can reach 
108 K; when this happens, these additional reactions must be in- 
cluded. They are I?N(p, y)!40, !?O(p,s)!9F, and F(p, y)!SNe. In 
each case, the product quickly decays via a positron emission: 140 — 
14N, 16F — 160. and 18Ne > 15F — 150, If any t80 is created, it 
will then combine with a proton and produce a *He nucleus via the 
reaction 150(p, a)?N ] 


A corollary to the above is that because the beta decay rates are 
so short, !?N, !?O, and !"F quickly reach equilibrium (in a matter 
of minutes). Thus, the time derivatives of (14.1.12), (14.1.15), and 
(14.1.18) can be set to zero, and these variables can be eliminated 
from the network. The calculation of the CNO cycle then reduces to 
solving six coupled differential equations: 
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Again, these equations can be solved numerically via Runge-Kutta 
integration. However, under most circumstances, the problem can 
be simplified enormously. Examine the table (again, from Clayton 
1968), which gives the reaction lifetimes as a function of burning 
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with p in g-cm ? and 7, in years. For the higher temperatures 
(which correspond to higher mass stars), the timescales listed are 
all much shorter than the timescale for the depletion of hydrogen, 
or for changes in the stellar temperature and density. Thus, while 
the CNO isotopes are achieving equilibrium, it is fair to assume that 
the hydrogen abundance and main sequence structure of the star are 


3 


constant. This is important, since it allows us to fix the values of 7 
in equations (14.1.26) - (14.1.31), and linearize the equations. 





To see this, let the subscripts 1 through 6 represent !?C, 13C, 14N, 
15N, 160. and 17O, and let y be the number the density of the species. 
In differential form for back-differencing, equations (14.1.20) through 
(14.1.25) are 
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These equations can now be re-written into a more useful form: 
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which, in matrix form is 
X 0 0 X 0 O0 V, Ui 
X X 0 0 0 0 y^ y2 
0 X X 0 0 X ys | luv 
0 0 X X 0 0 yl | Ya 
0 0 0 X X 0 Ys Ys 
0 0 0 0 X X Uo Yg 


where the X’s represent the terms of (14.1.32) - (14.1.37). The up- 
dated CNO abundances can therefore be found via a simple 6 x 6 
matrix inversion. 





The time it takes the CNO cycle to reach equilibrium is a critical 
function of temperature. However, two properties of the chain should 
be noted. First, the CNO equations really consist of two separate 
cycles: a CN cycle and an ON cycle. If you examine the lifetimes of 
the CN reactions, you can see that this cycle proceeds fairly rapidly, 
with the timescale for IN fusion (r ~ 10° years at Tg ~ 20 K) 
being the limiting factor. Since the branching ratio out of this cycle 
is very small (y ~ 4 x 1074), ~ 10? cycles are required before a 
significant number of CN nuclei are lost to the ON chain. Meanwhile, 
the timescale for ON nuclei to filter into the CN cycle is similarly long, 
T ~ 108 years (the !6O reaction timescale for Tg ~ 20 K). As a result, 
the CN cycle falls into quasi-equilibrium fairly quickly, and is largely 
decoupled from the slower ON cycle. 











The second thing to note is that in high mass stars (i.e., stars with 
high central temperatures), the reaction timescales are short, and 
the cycle achieves equilibrium fairly quickly. Thus, for purposes of 
estimating energy generation, the equilibrium abundances of the el- 
ements can (most times) be used. On the other hand, in low mass 
stars, where CNO equilibrium may never be achieved, the proton- 
proton chain dominates energy production. Thus, small errors in the 
CNO reaction are acceptable. 








Of course, the operation of the CNO cycle means that the differential 
equations for 'H and ^He must be modified to include the additional 
destruction and creation mechanisms. 


Table 5-3 Dependence of log (rp 4/100) on temperature t 


Temper- Reaction $ 
ature, | ———————————————————————————— 
T,  CU(p,y)N' C'(p,y)N'* Nis(p5)O? N'(pa)CH — 105 Op, y) FI OU(p,q)NH 
5 16.32 15.73 19.79 15.53 4.649 22.95 21.92 
6 14.32 13.73 17.57 13.29 4.598 20.51 20.02 
7 12.72 12.13 15.79 11,50 4.551 18.56 . 18.26 
8 11.41 10.81 14.32 10.03 4,508 16.95 16,50 
9 10.29 9.69 13.08 8.78 4,468 15.59 15.10 
10 9.33 8.73 12.02 7.70 4.431 14.42 14.05 
11 8.50 7.90 11.09 6.76 4.396 13.39 13.15 
12 7.75 7.15 10.26 5.93 4.363 12.49 12.38 
13 7.09 6,49 9.52 5.18 4.332 11.68 11.68 
14 6.49 5.89 8.86 4.51 4.303 10.95 11.02 
15 5.95 5.35 8.26 3.90 4.275 10.29 10.32 
16 5.45 4.85 71 3.34 4.248 9.68 9.55 
17 5.00 4.39 7.20 2.83 4.223 9.13 8.70 
18 4.58 3.97 6.73 2.35 4.198 8.61 1.86 
19 4.18 3.58 6.30 1.91 4.175 8.14 7.01 
20 3.82 3.21 5.89 1.50 4.1532 7.69 6.18 
22 3.16 2.55 5.16 0.75 4.110 6.89 4.78 
24 2.07 1.97 4.51 0.09 4.071 6.18 3.63 
25 2.30 1.70 4.21 —0.21 4.052 5.85 3.10 
26 2.05 1.44 3.93 —0.50 4.034 5.54 2.62 
28 1.58 0.97 3.41 —1.03 4.000 4.97 1.75 
30 1.15 0.54 2.93 —1,61 3.967 4.45 1.05 
35 0.23 —0.38 1.91 —2.55 3.803 3.33 —0.42 
40 —0.53 —]1.14 1.07 —8.42 3.829 2.41 —1.50 
45 —1.18 —]1.7$ 0.36 —4.14 3.771 1.64 —2.33 
50 —1.73 —2.33 —0.25 —4.77 3.719 0.97 —2.99 
55 —2.21 —2.82 —0.78 —3.32 3.073 0.39 —3.53 
60 —2.64 —3.24 —1.25 —5.81 3.630 —0.12 —3.97 
65 —3.02 —3.63 —1.67 —6.24 3.590 —0.58 —4.33 
70 —3.37 —3.97 —2.05 —6.63 3.554 —0.99 —4.65 
75 —3.68 —4.28 —2.39 —6.99 3.521 —1.37 —4.91 
80 —3.97 —4.57 —2.71 —1.32 3.489 —1.71 —5.14 
85 —4.23 —4.83 —2.99 —7.02 3.460 —2.02 —5.35 
90 —4.48 —5.08 — 3.26 —7.90 3.433 —2.31 —5.52 
95 ~4.70 —5.30 —3.51 -—-8.15 3.407 —2. 58 —5.68 
100 —4.91 —5.51 —3.74 —8.39 3.383 —2.83 —5.82 


OE UOUÁTÓ- 
1 Adapted from G. R. Caughlan and W. A. Fowler, Astrophys. J., 186:453 (1962). By permis- 
sion of The University of Chicago Press, Copyright 1962 by The University of Chicago. 


$ The lifetimes against protons are expressed in years, and the density p is in grams per cubic 
centimeter. ; 


Energy Generation in the CNO Cycle 


The energy produced by the CNO cycle is simply the sum from the 
individual reactions in the chain. However, we can approximate the 
energy generation in equilibrium by neglecting the ON chain, and 
considering that the “cycle time” for CN is dominated by the slowest 
reactions. Thus, we can approximate the chain with three steps: 


*C(p, y)'°N(B7, v) °C 


I3O(p, y) N 
V N(p, 4)? O(8* , v) N(p,* He)'?C 


The first reaction produces 1.944 + 2.221 — 0.710 = 3.455 MeV per 
reaction; the second generates 7.550 MeV per reaction, and the third 
14.019 MeV, per reaction. The total energy for the chain in equilib- 
rium is therefore 


N(?2C NGC N(14N 
712 T13 T14 
x 1079 ergs s^! cm~’ (14.2.1) 


or, since the abundance of the individual isotope is simply propor- 
tional to its lifetime. 
N(CN 
pec © 40.08 x 10 9 eC ae ergs s ! cm? (14.2.2) 
T12 + T13 + T14 

Now recall that these are non-resonant reactions. The nuclear life- 
times will depend on the reaction rates, A, which has a temperature 
dependence of 
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Note that this dependence is proportional to Z?/?. Thus, we can 
expect that the CNO cycle's energy generation will be ~ 5 times 
more sensitive to temperature than the proton-proton. 


A reasonable approximation to the CNO cycle's rate of energy pro- 
duction in equilibrium is 


cono = 8.67 x 10?” gia Xono X p Tg ^ e715228/T7^. ergs g-1 cm~ 
with 
9141 = 1 + 0.0027 T$? — 0.00778 T2/? — 0.000149 T; 


As Tg varies from 10 to 50, v goes from 23 to 13. This compares to 
v ~ A for the proton-proton chain. 
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Numerical Note: Forward and Back Differencing Techniques 


There are two ways to numerically solve a system of differential equa- 
tions. The first method is through “forward differencing,” or “ex- 
plicit” numerical integration. With this technique, values for inte- 
grating step n + 1 are estimated directly from the known values of 
the previous step. In the alternative “back differencing” or “implicit” 
integration method, the values at step n + 1 are estimated using the 
n + 1 values themselves, via numerical inversion. 





For example, consider the simple differential equation 


dX 

— =a-—bX 

di 
which has the analytical solution 

A m z Tpk 


In the forward differencing method, this equation is 


Anpi UU Xn 


Le 
At B 


which translates to 
Xn+1 = aAt+ Xn (1 — bAt) (4.2.1) 


while in the back differencing method, the equation is 


Apai m A 
— RAE. 08 EB 
P At X 
oaa (4.2.2) 


1+ bAt 


Now note that (4.2.1) and (4.2.2) behave very differently as the time 
step becomes large, i.e., At — oo. With forward differencing, the 
numerical solution diverges; with back differencing, X441 — a/b, 
which is the analytical result. This is characteristic of the two ap- 
proaches. In general, forward differencing solutions are unstable, and 
can only be used with very small time steps (i.e., for the above ex- 
ample, bAt < 1). Back differencing methods are much more reliable, 
and should be used whenever possible. 





